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CLAIMS 




1. A method of calculating the fast Fourier 
transform or the inverse fast Fourier transform of a 
5 digital signal defined by a series of N real starting 
samples\(n) , with N a power of two and n e [0..N-1], 
comprising^ successive transformation steps (2) for 
transforming^ input samples into output samples, all the 
transformation, steps being performed by means of a 
10 single set of\ butterflies with several inputs and 
several outputs\ the operating mode of which is 
modified selectively in each transformation step, the 

\ 

input and output samples of each transformation step 
being stored in a storage memory, a series of N output 
15 samples y(n) representative of the fast Fourier 
transform or the inverse, fast Fourier transform of the 



output samples x(n) beting provided in the last 

transformation step, 

characterized in that\ output samples y(n) are 

\ 

20 real, 

and in that the output \samples of a butterfly 

replace the corresponding input^ samples of the same 

rank in the storage memory, so that, if the starting 

samples x(n) processed in the f irst^transf ormation step 

25 are classified in bit-reversed orders of their index n, 

\ 

output samples y(n) are provided in the last 

\ 

transformation step in ascending order of index n, 
these output samples being defined by\the following 
relations : 

30 \ 



i^O) - Re[X(0) ] 

y(r^ = Re[X( (n+l)/2) ] for n being odd and 

different from N-1 
y(n) - x Im[X(n/2)] for n being even and 

different from 0 

y (N-1) -t Re [X (N/2) ] 
where samples\x(n) , with n e [0..N-1], designate the 
complex samples\ of the series corresponding to the fast 
or inverse fast; Fourier transform of the starting 

\ 

sample series x (nO . 

2, A method \ of calculating the fast Fourier 
transform or the inverse fast Fourier transform of a 
digital signal defined by a series of N complex samples 
X(n) conjugated by pairs represented by a series of N 
real starting samples \y (n) , with N power of two and n 

g [0..N-1], the starting samples y(n) being defined as. 

\ 

follows : 

y (0) = Re[X(0) ] 

y(n) = Re [X ( (n+1) /2) ] \ for n being odd and 

lifferent from N-1 
y(n) = Im[X(n/2)] for n being even and 

different from 0 

y(N-l) - Re[X(N/2) ] 




this calculation method comprising successive 
transformation steps for transforming input samples 
into output samples, a series of\sT output samples x(n) 
representative of this fast or Inverse fast Fourier 
transform being provided in the last transformation 
step, all the transformation steps being performed by 
means of a single set of butterflies with several 
inputs and several outputs, the operating mode of which 
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modified selectively in each transformation step, 
they input and output samples of each transformation 
step\being stored in a storage memory, 

characterized in that output samples x(n) are 

i \ 
real, \ 

\ 

and* in that the output samples of a butterfly 

\ 

replace tShe corresponding input samples of the same 

\ 

rank in the storage memory, so that, if the starting 

samples y (nk processed in the first transformation step 

\ 

are classifi-ed in ascending order of index n, the 

output samples x (n) are provided in the last 

\ 

transformation ^step in bit-reversed order of index n. 

\ 

3. The calculation method according to claim 1 or 
\ 

2, characterized \n that, in each transformation step, 
each butterfly transforms input sample pairs, the ranks 
of the input samples^ of the same pair within the series 



of input samples ofe said transformation step being 
symmetrical with respect to a center between the end 
rank values of the input samples transformed by said 
butterfly. \^ 

4. The calculation method according to claim 3, 

\ 

characterized in that it comprises u~l transformation 

\ 

steps E p with u=log 2 (N) and p e [0..p-2] . 

5. The calculation methcyd according to claim 4, in 



turn dependent on claim 3, iA turn dependent on claim 
1, characterized in further comprising: 

- a preliminary step of modifying the sequence of 



the starting samples x(n) ranked^in ascending order of 
index n and showing them in bit-reversed order of index 
n in the first transformation stepA and 
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- a final step of processing the series of output 

samples y(n) and providing a series of N complex 

conjugated samples X(n) corresponding to the fast or 

the inverse fast Fourier transform of the series of 

5 starting ^samples x (n) . 

6. The calculation method of claim 4, in turn 
\ 

dependent on claim 3, in turn dependent on claim 1, or 

according to\ claim 5, characterized in that, in each 

transformations step E p/ butterflies are distributed 

10 among N/2 P+2 calculation blocks, 

\ 

in that eath calculation block has a peripheral 



butterfly and/or v p -l internal butterflies, 

in that the\^peripheral butterfly of the rank a 
calculation block in transformation step Ep transforms 



15 the input samples o V f rank 2 p+2 ct, 2 p+2 a+2 p+1 -l , 2 p+2 a+2 p+1 , 

\ 

2 p+2 a+2 p+2 -l into output samples of the same rank, 

\ 

and in that the internal rank x butterfly of the 

\ 

rank a calculation block in transformation step Ep 



transforms the input ^samples of rank 2 p+2 a+2x+l, 
20 2 p+2 a+2x+2 , 2 p+2 a+2 p+2 -2x-3 , \^ 2 p+2 a+2 p+2 -2x-2 into output 
samples of the same rank, with (3>1 . 

7. The calculation metood according to claim 6, 
characterized in that each ^butterfly is assigned a 
coefficient W s , whereon the \calculation inside the 
25 butterfly is based, said coefficient being equal to e" 
j(2«s/N> W j_ th s G [0.. N/4-1] for a fast Fourier transform 
and is equal to e j(27ls/N) with s \f [0.. N/4-1] for an 



inverse fast Fourier transform. 



\ 



8. Calculation method according to claim 7, 
30 characterized in that the internal ra-nk x butterfly of 
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rank a calculation block in transformation step Ep 
assigned coefficient W s with 5 = (x+1) . (N/2 P+2 ) . 

9. The calculation method according to claim 8, 
characterized in that the butterflies for implementing 
the transformation steps are all of the same type and 
have \ 

four inputs for receiving input samples and 



four output>s for providing output samples, 

fourV additional inputs, respectively primary 
mode, secondary mode, permutation, and coefficient 



inputs, 
in 



order 




selectively apply different 



30 



transformation operations to the input samples, each 
operation being determined by the values assigned to 
the primary mode, secondary mode, permutation signals, 
and a coefficient admitted on said corresponding 
additional inputs . 

10. The calculation Method according to claim 9, 
characterized in that, for\each butterfly, the primary 
mode signal is 0 for a peripheral butterfly and 1 for 
an internal butterfly, \ 

in that the permutation saSgnal is 0 for the even 
rank calculation blocks, including rank 0, and 1 for 
the other ones. 

11. The calculation method according to claim 10, 
characterized in that, in transf ormauaon step Ep, each 
calculation block comprises one peripheral butterfly 
and 2 P -1 internal butterflies 

12. The calculation method according to claim 11, 
characterized in that the secondary mode sVgnal is 1 if 
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X 



e2 
e3 



and the 



i^e peripheral butterfly is used for the last 
transformation step, and otherwise 0. 

13. The calculation method according to claim 12, 
characterized in that, for four input samples el, e2, 
e3, ancke4, and for a complex coefficient W s =A+j*B, the 
butterf ly^delivers the following output samples si, s2, 
s3, and S ^S^ 

1) if the primary mode and secondary mode signals 

are 0: si =^el + e2 
\ 

s2 - el - 
s3 - e4 \ 
s4 = e3 

2) if the primary mode signal is 
secondary mode signal is 1 

si = el + e2 + \^3 + e4 
s2 = el - e2 
s3 = e4 - e3 

s4 = (el + e2) - (e\ + e4) 

\ 

3) if the primary mod^ signal is 1 and the 
permutation signal is 0: 

si = el + A.e3 - B.e4 
s2 = e2 + B.e3 + A.e4 
s3 = el - A.e3 + B.e4 
s4 = -e2 + B.e3 + A.e4 

4) if the primary mode signal is 1 and the 
permutation signal is 1 : 

si = el - A,e3 + B.e4 
s2 = -e2 + B.e3 + A.e4 
s3 = el - A.e3 - B.e4 
s4 = e2 + B.e3 + A.e4 
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14. The calculation method according to claim 10, 
characterized in that, in transformation step E p , each 
calculation block comprises: 

2 P -1 internal butterflies and a peripheral 
butterf ly\f or the even values of index p as well as for 
the last transformation step if p. is even, and 
2 p -u internal butterflies, otherwise. 

15. The \calculation method according to claim 13, 
characterized\in that the secondary mode signal is 1 if 

10 the peripheral butterfly is used for the last 

transformation step with y. being odd, and otherwise 0, 

\ 

16. The calculation method according to claim 15, 
characterized in that, for four input samples el, e2, 

\ <3 

e3, and e4, and fon a complex coefficient W=A+j.B, the 
15 butterfly delivers rhe following output samples si, s2, 
s3, and s4 

1) if primary mod^, secondary mode and permutation 
signals are 0; 

si = el + e2 + e% + e4 
20 s2 = el - e2 

s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 

2) if the primary moMe signal is 0 and the 
secondary mode signal is 1 

25 si = el + e4 

s2 = e2 
s3 = e3 
s4 = el - e4 

3) if the primary mode s\gnal is 0 and the 
30 permutation signal is 1: 

si = (e3 + e4) - (el + e2) 
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>2 = el - e2 
= e4 - e3 

= el + e2 + e3 + e4 

4) if \ the primary mode signal is 1 
permutation \signal is 0: 

si =\el + A.e3 - B.e4 
s2 = e2 + B.e3 + A.e4 
s3 = ei - A.e3 + B.e4 
s4 = -e2 + B.e3 + A.e4 

5) if the ^primary mode signal is 1 
permutation signal is 1 : 

si = el - A.e3 + B.e4 



and the 



and the 



s2 = -e2 + B.e3 + A.e4 
s3 = el + A.^e3 - B.e4 
s4 = e2 + B.eS + A.e4 

17. The calculation method according to claim 10, 
characterized in that, \ in transformation step E p , each 
calculation block comprises: 

2 P -1 internal butterflies and a peripheral 
butterfly for the even valoies of index p, and 
2 P -1 internal butterflies, otherwise. 

18. The calculation mekhod according to claim 17, 

\ 

characterized in that the secondary mode signal is 1 if 
the peripheral butterfly iv^ used for the first 
transformation step with \i being even, and otherwise 0. 

19. The calculation method\ according to claim 18, 
characterized in that, for four\ input samples el, e2, 
e3, and e4, and for a complex coefficient W s =A+j ,B, the 
butterfly delivers the following output samples si, s2, 



30 



s3, and s4 
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' x 1) vif the primary mode signal is 0 and the 
secondary mode signal is 1: 
ll = el + e2 



s2 = el - e2 
s3\= e4 - e3 
s 4 V= e3 + e 4 

2) if primary mode, secondary mode and permutation 
signals are 0: 

si = el + e2 + e3 + e4 
s2 = el - e2 
s3 = e4\- e3 
s4 = (eA + e2) - (e3 + e4) 

3) if the primary mode and secondary mode signals 
are 0 and the permutation signal is 1: 

si = (e3 +\e4) - (el + e2) 
s2 = el - 
s3 = e4 - eC 
s4 = el + e2\+ e3 + e4 

4) if the primary mode signal is 1 and the 
permutation signal is\0: 

si = el + A.e3\- B.e4 
s2 = e2 + B.e3 \ A.e4 
s3 = el - A.e3 +\B.e4 
s4 = -e2 + B*e3 +\A.e4 

5) if the primary \mode signal is 1 and the 
permutation signal is 1 

si = el - A.e3 + B.©4 
s2 = -e2 + B.e3 + A.\e4 
s3 = el + A.e3 - B.e< 
s4 = e2 + B.e3 + A.e4 
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yiL^ ^0 . The calculation method according to claim 8, 
characterized in that the butterflies for implementing 
the transformation steps are all of the same type and 
have 

5 - fomr inputs for receiving input samples and 

four outputs for providing output samples, 

- four\ additional inputs, respectively primary 
mode, secondary mode, permutation, and coefficient 
inputs, 

10 in order\ to selectively apply different 

transformation operations to the input samples, each 
operation being determined by the values assigned to 
the primary mode, \econdary mode, permutation signals, 
and a coefficient \ admitted on said corresponding 

15 additional inputs, 

and in that the f^inal step furthermore performs an 

addition and subtraction between the first and the last 

\ 

output sample provided in the last transformation step. 

21. The calculation method according to claim 20, 
20 characterized in that, in \transf ormation step E p/ each 

calculation block comprises one peripheral butterfly 
and 2 P -1 internal butterflies-. 

22. The calculation method according to claim 21, 

\ 

characterized in that, for four input samples el, e2, 

25 e3, and e4, and for a complex coefficient W s =A+j.B, the 

\ 

butterfly delivers the following Yputput samples si, s2, 
s3, and s4 \ 

1) if the primary mode signal is 0: 
si = el + e2 
30 s2 = el - e2 

s3 = e4 - e3 
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4 = e3 + e4 

2) \f the primary mode signal is 1 and the 
permutatibn signal is 0: 

sl\= el + A.e3 - B.e4 
s2 >f e2 + B.e3 + A.e4 
s3 A el - A.e3 + B.e4 
s4 = \-e2 + B.e3 + A.e4 

3) if th\ primary signal is 1 and the permutation 
signal is 1 

10 si = el\- A.e3 + B.e4 

s2 = -e2\ + B.e3 + A.e4 
s3 = el +\A.e3 -B.e4 
s4 = e2 + &.e3 + A.e4 
23. The calculation method according to claim 9 or 
15 20, characterized in\ that the first and second binary 
addresses of \i bits a\re generated for each butterfly, 
each binary address corresponding to the rank of an 
input sample of said butterfly and the second binary 
address being greater thcNj. the first binary address. 
20 24. The calculation n\ethod according to claim 23, 

characterized in that sai\i first and second binary 
addresses are consecutive arid an internal butterfly is 
involved. 

25. The calculation method according to claim 23 
25 or 24, characterized in that, i>f a peripheral butterfly 

is involved, the p+2 low-ord&r bits of the first 
address are equal to 0, and the\p+2 low-order bits of 
the second address form a number aqual to 2 P+1 -1. 

26. The calculation method according to claim 24 
30 or 25, characterized in that the \address of the two 

other samples to be applied to tihe inputs of the 
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itterfly, be they peripheral or internal, are obtained 
byXinverting the (p+2) low-order bits of said first and 
secqnd produced addresses . 

\27. The calculation method according to claim 26, 
characterized in that even-numbered address samples and 

\ , ■ 

odd-numbered address samples are stored in two separate 
memories^ 

28. The calculation method according to claim 25, 

\ 

characterized in that the value of the parameter s of 
\ 

the coefficient W assigned to an internal butterfly m 

transf ormatiorX step E p is coded by ]a-2 bits, and is: 

if p+l^yi-2, the number formed by the p+1 
\ 

low-order bits of the second binary address produced 

for said internal\butterf ly, 

if p+l<vi-2v the number formed by the p+1 

low-order bits of \he second binary address produced 

\ 

for said internal butterfly, followed by p-p-3 zero 



bits at the end of the ^number, 

if p+l>y-2, the\^n umber formed by the p+1 low- 
order bits of the second binary address produced for 
said internal butterfly, \ minus its p-p-1 low-order 
bits . 

29, The calculation method according to claim 4, 
in turn dependent on claim V3, in turn dependent on 



claim 2, characterized in that in each transformation 

step E p , the butterflies are^ distributed among 2 P 

calculation blocks, 

in that each calculation \block comprises one 

peripheral butterfly and N/2 P * 2 -1 internal butterflies, 

\ 

in that the peripheral butterfly of the rank a 
calculation block in transformation \step Ep transforms 
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khe input samples of rank 2 M " p a, 2 M ~ p a+2 M_p " 1 -l , 2 p_p a+2 p - p " 

x , \2 p " p a+2 p " p -l into output samples of the same rank, 

and in that the internal rank x butterfly of the 

rank a calculation block in transformation step En 
\ 

transforms the input samples of rank 2 U p a+2x+l, 2 M 
p a+2x+2, \^~ p a+2 p ~ p_1 -2x-3, 2 M ~ p a+2 M_p -2x-2 into output 
samples of the same rank. 

30. The^calculation method according to claim 29, 

characterized \ln further comprising a final step of 

modifying the sequence of the output samples provided 

\ 

in the last transformation step and classifying them m 
ascending order of index n. 

31. The calculation method according to claim 29 

\ 

or 30, characterized in that each butterfly is assigned 
a coefficient W s , whereon the calculation inside the 
butterfly is based, saidNcoef f icient being equal to e" 
j(2K S /N) with s € [0. .N/4-1] Tor a fast Fourier transform 

and is equal to e j(2ns/N) wxth s e [0.. N/4-1] for an 

\ 

inverse fast Fourier transform. 

\ 

32. Calculation method according to claim 31, 

\ 

characterized in that the interrfal rank x butterfly of 
the rank a calculation block in transformation step Ep 
is assigned coefficient W 5 with 5 = - 

33. The calculation method according to claim 32, 
characterized in that the butterf liesVf or implementing 
the transformation steps are all of the same type and 
have 

- four inputs for receiving input \ samples and 
four outputs for providing output samples, 



o 9 s s ± s: 7 1 e ,„ in, en :i 7' o e 



49 




10 



15 



20 



25 



30 



four additional inputs, respectively primary 
secondary mode, permutation, and coefficient 



order 



to 



selectively apply different 



transformation operations to the input samples, each 
\ 

operation being determined by the values assigned to 
the primary\mode, secondary mode, permutation signals, 
and a coefficient admitted on said corresponding 
additional inputs . 

34. The ca\culation method according to claim 33, 
characterized in Vhat, for each butterfly, the primary 
mode signal is 0 for a peripheral butterfly and 1 for 
an internal butterfly, 

in that the permutation signal is 0 for the even 
rank calculation block v s, including rank 0, and 1 for 
the odd values. 

35. The calculation \method according to claim 31 

or 34, characterized in that the secondary mode signal 

\ 

is 1 if the butterfly, be it^peripheral or internal, is 
used for the first transformation step, and otherwise 



0. 



\ 



36. The calculation methods according to claim 35, 
characterized in that, for four\ input samples el, e2, 
e3, and e4, and for a complex coefficient W s =A+j.B, the 
butterfly delivers the following output samples si, s2, 



s3, and s4 

1) if the primary mode and sec 
are 0 : 

si = (el + e2)/2 
s2 = (el - e2) /2 
s3 = (e4 - e3)/2 



idary mode signals 
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s4 = (e3 + e4) /2 

if the primary mode signal is 0 and the 

secondary mode signal is 1: 

si = [ (el + e4) /2-e2] /2 

5 s 2 = [ (el+e4) /2-e2] /2 

s^v= [e3- (el-e4) /2] /2 

s4 = [e3+ (el + e4) /2] /2 
\ 

3) if the primary mode signal is 1 and the 

\ 

permutation signal is 0: 
10 si = (W + e3)/2 

s2 = (e2+e4) /2 

s3 = [(el-e3).A - (e2+e4).B)/2 

s4 = [-(el-e3).B + (e2+e4).A]/2 
\ 

4) if the prvimary mode signal is 1 and the 
15 permutation signal rs 1: 

si = [ (el-e3)\A - (e2+e4).B]/2 
s2 = [-(el-e3)\B + (e2+e4).A]/2 
s3 = (el+e3) 12 
s4 = (e2-e4)/2 

20 37. The calculation method according to claim 33, 

characterized in that the first and second binary 
addresses of p bits are generated for each butterfly, 
each binary address corresponding to the rank of an 
input sample of said butterf\y and the second binary 

25 address being greater than the first binary address. 

38. The calculation method according to claim 37, 
characterized in that said firs\t and second binary 
addresses are consecutive and an internal butterfly is 
involved. 



\ 



30 39. The calculation method according to claim 37 

\ 

or 38, characterized in that, if a peripheral butterfly 



involved, the p-p low-order bits of the first 

adaress are equal to 0, and the p-p low-order bits of 

the second address form a number equal to N/2 P+1 -1. 

4(X The calculation method according to claim 38 
\ 

or 39, characterized in that the address of the two 

other samples to be applied to the inputs of the 
\ 

butterfly are obtained by inverting the p-p low-order 
bits of both produced addresses. 



41. The calculation method according to claim 40, 

\ 

characterized in tl^at even-numbered address samples and 
odd-numbered address^ samples are stored in two separate 
memories . 

42. The calculation method according to claim 41, 
characterized in that the value of the parameter s of 
the coefficient ttr assignedv^to an internal butterfly in 

transformation step E p is coded by p-2 bits, and is: 

\ 

- if p-p-l=p-2, the number formed by the p-p-1 

low-order bits of the second address produced for said 

internal butterfly, 

if p-p-Kp-2, the number ^formed by the p-p-1 

\ 

low-order bits of the second address, produced for said 
internal butterfly, followed by p-1 \zero bits at the 
end of the number, 

- if p-p-l>p-2, the number formed^ by the p-p-1 
low-order bits of the second address produced for said 
internal butterfly, minus its p+1 low-order\bits . 



